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Abstract The creation of point defects in the crystal lattices of various semicon¬ 
ductors by subthreshold events has been reported on by a number of groups. These 
observations have been made in great detail using sensitive electrical techniques but 
there is still much that needs to be clarihed. Experiments using Ge and Si were 
performed that demonstrate that energetic particles, the products of collisions in the 
electron beam, were responsible for the majority of electron-beam deposition (EBD) 
induced defects in a two-step energy transfer process. Lowering the number of col¬ 
lisions of these energetic particles with the semiconductor during metal deposition 
was accomplished using a combination of static shields and superior vacuum result¬ 
ing in devices with defect concentrations lower than 10^* cm^^, the measurement 
limit of our deep level transient spectroscopy (DLTS) system. High energy electrons 
and photons that samples are typically exposed to were not influenced by the shields 
as most of these particles originate at the metal target thus eliminating these parti¬ 
cles as possible damage causing agents. It remains unclear how packets of energy 
that can sometimes be as small of 2 eV travel up to a /rm into the material while 
still retaining enough energy, that is, in the order of 1 eV, to cause changes in the 
crystal. The manipulation of this defect causing phenomenon may hold the key to 
developing defect free material for future applications. 
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1 Introduction 

Process induced defect creation in semiconductors is of paramount importance as 
device performance is influenced, adversely or beneficially, by these defects ll2^ . 
Semiconducting materials offer the ideal platform for studies into point defects with 
energy levels in the bandgap as ultra-pure material is readily available and can then 
be investigated using techniques like deep level transient spectroscopy (DLTS) 
to measure the energy level of the defect, also known as the defect enthalpy, de¬ 
fect concentration and apparent capture cross-section. Additionally, Laplace DLTS 
resolves two or more defect levels that present as a single broad peak in the con¬ 
ventional DLTS spectrum m further clarifying complex observations. While these 
techniques are unable to provide a physical description of a defect, they are sensitive 
to defect concentrations as low as 10^' cm^^, in our experiment. Semiconductors 
are technical materials that now enable us to directly measure the effects of radiation 
on structured systems. 

A sample is typically exposed to 10 keV electrons during electron-beam deposi¬ 
tion (EBD) although sources with higher acceleration do exist, none of them exceed 
60 keV. All electron beam (EB) heated sources rely on energy transfer from incident 
electrons to thermally evaporate any one of a large variety of solid targets. The mod¬ 
ern electron gun (EG or E-gun) that was introduced in the early 1960s, remaining 
virtually unchanged since then, has found application in metallization on semicon¬ 
ductors, optics II 22 I and in industrial processes like the deposition of corrosion pro¬ 
tective coatings on strip metal ll35l . A detailed description of the EBD source and 
the power supplies that control it ll22l will not be repeated herein, however a brief 
description will follow to describe the source used for this investigation. 

The EBD source consists of three components: the electron emitter, magnetic 
lens and water-cooled cavity or hearth. The emitter is strategically located out of 
line-of-sight of the evaporant and the electron beam follows a circular path curved 
by the magnetic lens through 270° so as to impinge on the centre of the hearth. This 
protects the emitter from becoming coated by the evaporant, thus lowering the risk 
of short circuits and also conveniently shields the substrate from energetic particles 
that may be accelerated by the high potential of the emitter. Three power supplies 
are required, first to heat the filament (tungsten coil) thus providing a source of elec¬ 
trons, secondly to accelerate these electrons and finally to power the electro-magnets 
of the lens to control the electron beam. 10 kV is the most common accelerating 
voltage at a current of up to 1.5 A and was the source used for this investigation. 
In modern systems most tetrode based high voltage power supplies have been re¬ 
placed with solid state equivalents that are well protected from short circuits due to 
arcing. Eor safety in operation, today’s electron guns have a magnetic lens that con¬ 
sists of a permanent magnet to direct the electron beam towards the hearth centre, 
as well as electro-magnets to focus and raster the beam. Modern magnet supplies no 
longer defocus the beam to cover a larger area of the evaporant but rather maintain 
a focused beam that is scanned over the target surface in a complex pattern at a fre¬ 
quency not exceeding 200 Hz. This arrangement ensures that the target material is 
evenly heated thus better utilised and should the magnet supply fail then the electron 
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beam remains focused on the centre of the hearth. During operation efficient water 
cooling is of paramount importance if the hearth is to remain inert so as to ensure 
the purity of the deposited film. 

The disadvantage of EBD is that it introduces defects in sensitive semiconductors 
ll9l lT7ll28]l . This damage has previously been attributed to an emission of soft x- 
rays or energetic electrons that are most probably reflected from the target ||22|. 
The magnetic field of the E-gun will cause the majority of reflected electrons to 
be captured by the shield placed over the permanent magnet and is a signihcant 
part of the design as approximately 30% of the beam energy is reflected. A small 
portion of the evaporant flux is ionised as it passes through the incident electron 
beam further complicating matters. Another source of energetic particles that has 
previously been neglected is those ions that are created in the electron beam path 
by collisions between electrons and residual gas atoms or molecules. Even for fast 
moving atoms like hydrogen the probability of collision while traversing a typical 
10 kV, 0.1 A electron beam is above one. Eurthermore, as an evaporation proceeds 
the vacuum pressure tends to increase with increasing outgassing due to heating of 
the vacuum chamber and the components in the chamber, resulting in the number 
of available particles that may undergo collisions increasing proportionately with an 
increase in pressure. 


2 Radiation enhancement through intermediate collisions 

To account for subthreshold electron damage a two-step process was suggested 0 
|T5][3Tl[33l where an intermediate light impurity atom, such as hydrogen, could pro¬ 
duce a displacement of a germanium atom. This process requires the electron to 
first strike the light atom that then strikes the germanium atom transferring almost 
three times more energy than a direct collision. The electron threshold energy for 
such a displacement was found to be 90 keV, assuming that 15 eV is required to dis¬ 
place a germanium atom from the lattice IfT^ . While this threshold is much higher 
than the typical available electron energy, defects observed in gold and copper were 
postulated to be due to ever present impurity atoms cni. Similarly, in germanium, 
light-atom impurities are the most probable subthreshold mechanism agent. Naber 
and James ll^ only considered atoms present in the crystal lattice, but using light 
atoms that are present in the vacuum to transfer energy to lattice atoms theoretically 
yields the same result. Erom conservation of momentum and energy, if we consider 
two particles denoted by the subscripts 1 and 2 then let m\ and m 2 be the masses, 
u\ and U 2 be the velocities before collision and vi and V 2 be the velocities after an 
elastic collision then; 


m\U\ +m2U2 = niivi +m2V2 


( 1 ) 


and 



1 2,1 2_1 2,1 


(2) 
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1 10 

Atomic mass (AMU) 


Fig. 1 Theoretical maximum energy transfer in an elastic collision of a 10 keV electron (black plot 
with filled circles) or between a 24 eV hydrogen atom (blue dash plot) and particles of increasing 
mass. Bars denote the energy variation dependent on the velocity of the second particle in a vac¬ 
uum, parallel to the direction of the impinging particle. Relativistic considerations were included 
but only accounted for a 0.9% increase in energy transferred. 



Fig. 2 Maximum energy transfer in an elastic collision between a particle (particle 1) and a sta¬ 
tionary Ge lattice atom (green plot) where particle 1 was initially accelerated by a collision with a 
10 keV electron (black solid plot). The dashed line represents the maximum energy that is trans¬ 
ferred irrespective of particle 1 mass. A maximum energy of 1.3 eV was transferred to Ge using H 
as the intermediate atom whereas directly, only 0.34 eV was transferred 
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For the simplest case of U 2 = 0 the maximum energy transferred to particle 2 is 
given by; 


K,, 


1 2 1 2 4mim2 Am\m2 

-m2V2 = 2'”‘“'(mi+m2)2 ^ ^\mx+m2Y 


(3) 


assuming a one dimensional case of an elastic collision where £, is the initial en¬ 
ergy of particle 1. This energy transfer between a 10 keV electron and particles of 
atomic mass from 1 to 75 is illustrated in Fig.[T](black solid plot with filled circles). 
The red bars denote the energy variation if the velocity of the second particle in 
vacuum is taken into account and including this consideration then the maximum 
energy transferred to a H atom is approximately 24 eV. Plotting the example of col¬ 
lisions between a 24 eV H atom and particles of atomic mass 1 to 75 illustrates that 
this knock-on process is capable of transferring the same (only for AMU = 1) or 
more energy than a direct collision with an electron. To evaluate this process for the 
specific case of Ge, Figure |2] plots the knock-on energy transfer between particles 
of various masses that were initially accelerated in a 10 keV electron collision and 
then collide with a stationary Ge atom. Collisions of the lightest particles with Ge 
result in the highest energy transfer, that is, at most, 1.3 eV. This is not sufficient 
to displace a Ge atom from the lattice but Chen et al na noted that defects were 
only produced in Ge grown in a H atmosphere thus it is likely that H in the crystal 
lattice played a role. The direct electron-Ge elastic collision process only resulted 
in 0.34 eV being transferred to a stationary Ge atom. 

It is also of interest to know how the transfer of energy in collisions depends on 
the angle. Supposing that a particle of mass mi and kinetic energy K\ experiences a 
collision with a particle of mass m 2 and this one exits the collision with energy K 2 
with angle 02 with respect to direction of the incident particle. Then, it is easy to 
demonstrate that the curve K 2 / max(/r2) with respect to 02 does not depends on the 
masses or on the energy of the incident particle as can be seen in Fig. [3 There is a 
significant interval of exit angles for which K 2 is close to the maximum. 


Fig. 3 Relative kinetic energy 
with respect to the maximum 
possible kinetic energy of a 
particle that has experienced 
a collision while at rest, as a 
function of the exit angle with 
respect to the incident particle 
direction. This curve neither 
depends on the masses nor on 
the initial energy. It can be 
seen that the there is a large 
interval of angles for which 
the final energy is close to the 



Angle 


maximum 
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3 EBD experimental details 


A Ge (111) wafer, bulk grown and doped with Sb to a concentration of 1 x 10*^ 
cm^^ was degreased in successive 5 minute ultrasonic baths of trichloroethylene, 
isopropanol and methanol before being etched for 1 minute in a solution of 5;1 
H20;H202 (30%). To create an ohmic contact, AuSb was deposited resistively on 
the wafer back surface and then annealed in an Ar ambient at 350°C to lower the 
contact resistance. Samples cut from this wafer were then degreased and etched 
again before EBD of Pt through a metal contact mask was carried out to yield eight 
Schottky barrier diodes (SBDs) with a diameter of 0.6 mm and 50 nm thick on each 
sample’s front surface. All SBD depositions were carried out using an electron beam 
with an accelerating voltage of lOkeV and beam current of approximately 100 mA. 
Current-voltage measurements were carried out on all diodes to verify their suitabil¬ 
ity for DLTS analysis. 

Conditions in the EBD chamber were varied during diode manufacture by a) 
not applying any counter measures, b) back-filling the chamber with forming gas 
(H 2 :N 2 , 15%:85%) to 10^^ mbar, c) back-filling with forming gas and placing one 
shield (Bi in Fig.|4|i to shield from direct particles, d) back-filling with forming gas 
and applying 2 shields (Bi and B 2 ) so that particles reflected off the chamber wall 
are also shielded for and e) superior vacuum with low H 2 concentration as well as 
both shields in place. The measures taken to ensure that the H 2 concentration was 


Fig. 4 EBD chamber layout 
detailing the po.sitioning of 
static shields (Bi and B 2 ) 
used to shield samples from 
energetic particles created 
in collisions with the high 
energy electrons of the elec¬ 
tronic beam (C) 
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maintained below 10 ^ mbar and the DLTS spectra obtained have been published 
previously in). 

To investigate the role of energetic particles arriving at the semiconductor sur¬ 
face during EBD, clean samples were exposed to the conditions of EBD without 
any evaporation taking place, termed electron beam exposure (EBE) herein, and 
thereafter Schottky barrier diodes (SBDs) were evaporated resistively onto the ir¬ 
radiated Ge. These samples were exposed for 50 minutes at 100 mA beam current 
as this was approximately the same amount of exposure that the Pt EBD diodes 
received. The same measurement procedure was followed as applied previously. 


4 Defects after EBD and their origin 


The defects introduced during EBD have been reported on before of which the 
E-center is the most prominent. This defect consists of a vacancy-dopant complex, 
the dopant in this case being Sb. A control sample manufactured using resistive 
evaporation RE), a technique known not to introduce defects in Ge, had no mea¬ 
surable defects in it. The peak heights of the DLTS spectra are indicative of defect 
concentration as 


Nt _ 2AC 


(4) 


where Nt is the deep level concentration. No is the concentration of shallow impuri¬ 
ties, AC is the DLTS peak height and C is the junction capacitance. The capacitance 
of all the devices manufactured was found to be approximately the same and thus 
spectra can be compared directly. 

The DLTS spectra in Eigs. |5] and |6] were all obtained from diodes prepared in 
the same EBD system. Eor spectrum a) a standard oil-filled rotary vane pump was 
used but for all the other spectra an oil-free pump was used as the fore-pump during 
deposition. To further improve the vacuum all crucibles were baked out in situ us¬ 
ing the electron gun. It is important that the pressure not increase drastically during 
evaporation although a change in vacuum pressure is inevitable as fixtures heat up 
during EBD and then outgas. Comparing spectrum a) with spectrum b) it is evident 
that the peak heights of all the defects that are present in both spectra are reduced 
by approximately 90% in spectrum b). A further reduction in peak heights can be 
observed in spectra c), d) and e) once shields were applied. Spectrum e) that repre¬ 
sents a diode prepared in a superior vacuum with two shields in place presents as a 
wavy plot, indicative of surface states, but sharp peaks that are evidence of defects 
with deep levels are conspicuously absent. Shields Bi and B 2 were only capable 
of blocking off energetic particles that were created when 10 keV beam electrons 
collide with residual gas atoms or molecules and not for electrons reflected off the 
evaporant surface. Also, it is expected that light ions will follow a curved trajectory 
around shield Bj while acted on by the magnetic field of the electron beam thus 
rendering the shield ineffective to some degree. 
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Fig. 5 DLTS spectra recorded after electron beam deposition of Pt Schottky barrier diodes un¬ 
der various vacuum conditions. For plot a) standard vacuum conditions apply and there were no 
shields. For samples b), c) and d) the chamber was back-filled with forming gas at a partial pres¬ 
sure of 10“^ mbar where b) no shields, c) one shield (direct particles) and d) two shields (also 
the reflected particles) were applied. Plot e) represents a diode evaporated in a superior vacuum 
with two shields in place. Reproduced with permission from Coelho, S.M.M., Auret, F.D., Janse 
van Rensburg, P.J., Nel, J.: Electrical characterization of defects introduced in n-Ge during elec¬ 
tron beam deposition or exposure. J. Appl. Phys. 114(17), 173,708 (2013). Copyright 2013, AIP 
Publishing LLC. 


The large difference in defect concentration between spectra a) and b) was sur¬ 
prising when one considers that the only difference in the conditions was that less 
hydrocarbon contamination was present during the manufacture of sample (b) and 
that the vacuum pressure was kept constant at 10^^ mbar by introducing forming 
gas into the chamber. Sample (a) was initially at a vacuum pressure of 10^® mbar 
when the deposition started but this pressure quickly increased to 10^^ mbar or 
more as the chamber heated up. The composition of the residual gas present during 
EBD appears to be the largest contributor to the high defect concentration in sam¬ 
ple (a) as well as a higher partial pressure near the electron gun where most of the 
outgassing occurs. Crucibles used during these evaporations may also have played 
a role as the standard carbon crucible that was used when preparing sample (a) had 
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Fig. 6 DLTS spectra of a) RE Pd diode - the control, b) Pd EBD diode, c) , d) and e) EBE diodes 
prepared by RE. Diode c) received 50 minutes of EBE followed by RE Pd. Diode d) was exposed to 
5x10 minutes EBE followed by RE Pd and diode e) received 5x10 minutes of EBE followed by 
RE Au. DLTS measurement conditions were as stated on the figure. Reproduced with permission 
from Coelho, S.M.M., Auret, F.D., Janse van Rensburg, P.J., Nel, J.: Electrical characterization of 
defects introduced in n-Ge during electron beam deposition or exposure. J. Appl. Phys. 114(17), 
173,708 (2013). Copyright 2013, AIP Publishing EEC 


a greater impact on the vacuum pressure than the Fabmate® crucible that was used 
for other samples. The complex nature of conditions present during EBD is evident 
in Figs|5]and|6]if one considers that the defect concentration increased slightly with 
the addition of a second shield, a measure designed to lower the defect concentra¬ 
tion. This small difference was however not enough to draw conclusions from but 
most important was that all the counter-measures together lowered the defect den¬ 
sity to a level that could no longer be measured. For a diode evaporated onto Si the 
same measures lowered the defect density so that, although the DFTS peaks were 
small, some defects could still be identified Q . 

The maximum energy that can be transferred by impinging atoms to the Ge lat¬ 
tice per collision can be read off graph |2] (thick green plot). Treating all collisions 
elastically is a reasonable simplification to obtain the maximum possible energy 
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that can be transferred whereas treating the electrons relativistically only served to 
increase the energy transferred by 0.9% and need not be taken into account. The 
maximum energy that can be transferred to Ge was found to be approximately 1.3 
eV via a light atom like H, with maximum energy transferred decreasing as the 
intermediate atom or particle increases in mass. This energy is not sufficient to dis¬ 
lodge a Ge atom from its position in the crystal lattice but may dislodge a light 
atom that has taken up a substitutional position in the lattice or modify an existing 
defect with an energy level too close to the band edges to be detected using DLTS. 
Vacancy-hydrogen complexes have been reported on previously in Ge HI and were 
observed using infrared spectroscopy HD- There is at present no certainty which of 
these complexes plays a role in defect formation during EBD as their concentration 
in bulk grown Ge is too low to be detected with infrared spectroscopy. 

The defects that were observed after electron beam exposure of Ge that numbered 
ten different defects, in total, have not been observed before with the exception of 
Eo .37 and Eq ss (E-center). The defect concentration of EBE induced defects was 
much lower than that measured after EBD for similar exposure times and this is 
evident if one compares plots b) and c) of Eig. |6] During EBD the semiconductor 
receives a measure of protection from impinging particles as it is exposed to radi¬ 
ation through an ever increasing metal film. No such layer is present during EBE 
thus it was expected that similar or more damage would be observed after the EBE 
process. One possibility for the great variety of different defects observed is that 
these defects are mostly due to atoms being implanted into the EBE treated Ge but 
this cannot explain the absence of the EBD induced defects. That the metal layer 
acts as a channel for energy to be transferred to the semiconductor is a possibility 
that will require further investigation. Samples exposed for 50 minutes in 10 minute 
increments interrupted with 50 minute periods to allow for cooling exhibited signif¬ 
icantly higher defect concentrations for all defects observed. The sample that was 
subjected to a continuous 50 minute EBE heated up 35°C more than the sample that 
was allowed to cool. Differences in defect concentrations may be due to annealing, 
in part, but cannot explain why all the EBE induced defects were equally affected. 
The other possibility is that sample heating interrupts the energy transfer process 
leading to less defects being introduced. Detailed annealing studies will be required 
to shed more light on this result. 


5 Intrinsic localized modes and defects 

Eor many years the paradigm of considering phonons as the entity transporting en¬ 
ergy in a solid has been overwhelming. Phonons as it is well known are obtained 
under the hypothesis of small lattice vibration that allows the linearization of the 
dynamical equations of the system or equivalently allows the use of the harmonic 
approximation for potentials. Perhaps one of the clearest example of success was 
Einstein solid theory where phonons were quantized in Ref. lISTl at the beginning of 
XXth century. Linear systems and phonons have been extremely successful not only 
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in the framework of quantum mechanics but also classical mechanics linear lattice 
theory has been very productive. Most of the theory of spectroscopy is based on the 
harmonic approximations and phonons. 


5.1 Limitations of harmonicity 

It is however based on several assumptions that are known to be convenient mathe¬ 
matical tools but not accurate representations of reality. First, it is well known that 
interatomic potentials are not harmonic, starting from the electrostatic interaction 
and continuing from Van der Waals forces described for example with Buckingham 
potentials V = Aexp(—E/kBT) —Bjr^. However, the harmonic approximation is 
quite convenient at temperatures of the order of room temperature and above, for 
which the average atomic displacement are not too large. The key word is aver¬ 
age, for average displacements or properties. At any temperature there is a small 
but finite probability that some displacements are large enough for the harmonic 
approximation to become invalid, but they will have a small effect in the average 
properties. However, even if considering only bulk properties, it is well known that 
the harmonic approximation is not sufficiently accurate as such a solid would not 
experience thermal expansion and would have an infinite thermal conductivity Q. 

There is a huge change, when interaction with radiation or swift particles is con¬ 
sidered. In this article, for example, we considered the possible interaction of very 
low energy particles such as lOkeV electrons or 24eV H atoms. Germanium atoms 
may acquire energies of 1 eV, forty times larger than the average thermal energy 
at room temperature. For the displacements involved we can be sure that nonlinear 
effects will take place. If the interatomic distances become small enough, potentials 
with a strong repulsive core such as Lennard-Jones or Ziegler-Biersack-Littmark 
(ZBL) 041II42II need to be introduced to provide a realistic description of the forces. 
If the energies are large enough they will produce defects in the solid by displac¬ 
ing atoms from their lattice positions, bringing about the formation of point defects 
like interstitials or vacancies. In this article and in this section we will focus our 
attention in energies that are not large enough to disrupt the lattice geometry, the 
so called subthreshold radiation regime. The threshold energy in Ge depends on the 
lattice direction, being 11.5eV and 19.5eV for the (111) and (100) directions, re¬ 
spectively 1251. Conventional knowledge supposes that the energy just disperses into 
phonons elevating locally the temperature of the sample creating a thermal spike but 
which soon would relax to thermal equilibrium with the rest of the crystal. 

Another shortcoming of the phonon description is that phonons are harmonic 
waves that extend over the whole space. This is a very useful mathematical hypoth¬ 
esis and it is justified because the extension of the phonons is much larger than the 
lattice unit. However, the impact of a lOkeV electron or a 24eV H atom on Ge is 
clearly a localized phenomenon, because the de Broglie wavelength is ^ 10^^ nm, 
smaller than an atom size. In the harmonic approximation the consequence of such 
an impact is a wave packet but because basically all media are dispersive it soon dis- 
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perses into phonons with different wavelengths and velocities and the localization 
is lost. 

However, if the nonharmonicity of the potentials is taken into account such an 
impact may produce what is called an intrinsic localized mode (ILM), also known 
as a breather, depending on the context 1^1^[^[^ . This is a localized wave 
packet that does not spread, that is, it behaves like a quasiparticle. As ILMs are not 
exact solutions they will eventually lose energy and disperse into phonons. How 
long they can live, how many of them are there and how important they are, are still 
open questions that are very much related and that we address here briefly. The key 
concept to understand breather existence is the fact that the frequency of vibration 
of nonlinear oscillators depends on the amplitude or energy of them, which does 
not happen in a linear oscillator. If the frequency of the oscillator increases with the 
amplitude, it is called a hard potential. This corresponds to a potential that grows 
faster with the distance to the equilibrium point than what the harmonic one does, 
while being equal at small distances. If the frequency of the oscillator decreases with 
energy, it is called a soft potential and it grows more slowly than the harmonic one 
does, as is illustrated in Fig. |7] The phonon spectrum of a solid is always bounded 
from above, may have gaps, and in some cases may also be bounded from below, 
in which it is called optical. If it is not bounded from below it is called acoustic. 
Vibrations with frequencies that are outside the phonon spectrum cannot propagate 
in the solid, bringing about localization of energy that does not spread. Figure [8] 
shows an example for a model of cations in a silicate layer which produces the 
optical spectrum El. 

Breathers are well described mathematical objects and are easy to produce in 
macro and meso systems. For example, a chain of magnetic pendulums is easy to 
construct and to experiment with liTTl . Their existence in a solid is a more difficult 
question for several reasons, to cite a few: a) the reality is quantum and not clas¬ 
sical; b) the potentials are simplifications of complex interactions; c) the lattice is 
not perfect; d) the lattice is disordered due to temperature. These subjects have been 
studied, theories of quantum breathers exist and molecular dynamics using increas¬ 
ingly realistic and complex potentials have been useful in creating ILMs of energies 
of the order of magnitude of eV O23ll24ll40l that propagate at finite temperatures. 
But more importantly, there is growing experimental evidence of long range local¬ 
ized transmission of energy. For example, it was observed ll^ that subsequent to 
the impact of an alpha particle on the surface of an insulator, there was transmission 
of energy in a localized way along close packed lines that was able to eject an atom 
at the surface of the crystal. For the material of interest in this article, germanium, 
it was shown that the impact of Ar atoms of 2-8 eV were able to anneal defects 
2 /rm below the surface GM- Annealing and ordering of voids in several crystals 
attributed to ILMs ifT^ is just another example. 
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Fig. 7 Left: harmonic potential (continuous line), soft potential (dots) and hard potential (dashed). 
Right: dependence of the frequency with respect of energy for an oscillator with different potentials 
(same code) 




E (KJ/mol) E (KJ/mol) 

Fig. 8 Examples of phonon spectra for atoms in a system with a substrate which produces an opti¬ 
cal phonon spectrum bounded from below. Left: the potentials are soft, and therefore, the energies 
of the ILMs diminished below the phonon band. Right: potentials are hard and the frequency of 
ILMs increase with energy. This was obtained for a model of cations in a silicate layer. Reproduced 
with permission from: Archilla, J.F.R., Cuevas, J., Romero, F.R.: Effect of breather existence on 
reconstructive transformations in mica muscovite. AIP Conf. Proc. 982(1), 788-791 (2008). Copy¬ 
right by American Institute of Physics 
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5.2 Ejfect of intrinsic localized modes 

One question is if ILM do exist in a solid, what will be the effect and on which prop¬ 
erties. If there are many of them they will probably interact between them and will 
be dispersed. The main effect will be an increase of the temperature of the system, 
similarly for a harmonic lattice. It seems that one of the most important effects of 
breathers could appear in connection with changes of structure, annealing, chemical 
reactions and similar processes. Generally speaking all processes for which a poten¬ 
tial barrier with some activation energy Eg has to be overcome and with a probability 
of happening proportional to exp{—Ea/kBT), that is, the constant rate of the process 
is given by an Arrhenius type equation: 

K=Aexp{-Ea/kBT). (5) 

This equation is extremely sensitive to changes in £« and it is also asymmetric, i.e., 
the increase in the rate K corresponding to a decrease of energy AE is much larger 
than the decrease in the rate corresponding to an increase of the same amount of 
energy. An easy calculation shows it. Suppose that there is some perturbation of the 
barrier AE during some time At and a perturbation AE during the same time, then, 
the mean rate k' during the time interval 2At would be: 

^_L(Ae-(Ea-AE)/kBT^^^^^-iE, + AE)/kBT^A 
2At V / 

^l_I^^AE/ksT ^^-AE/ksTy^-EjkBT 

The amplification factor is / = cosh{AE/khT) and can usually be approximated by 
I — j exp{AE/ksT). It does not depend on the height of the barrier Ea but only on 
the ratio of the barrier variation AE and the thermal energy of the lattice. It can be 
seen in Fig. |9] An elaborate and rigorous theory is developed in Refs. O17lll8ll20l . 
but the conclusions are the same. Therefore, ILMs of small energy, both mobile or 
stationary can produce a huge effect. Even more if we consider that their energy is 
localized and not extended as for phonons. 

An example of this phenomenon in a silicate is described in Ref. E). In an ex¬ 
periment, reconstructive transformation of the mica muscovite into lutetium disili¬ 
cate was observed to occur several orders of magnitude faster than expected due to 
the nature of the bonds that have to be broken. The explanation is based on a fact 
observed in numerical simulations: that breathers with larger energy live longer, 
therefore a temporary fluctuation that produces an accumulation of vibrational en¬ 
ergy and creates an ILM is not immediately destroyed ll34l . The more energetic the 
ILM the more unlikely, but also the longer the lifetime. Eventually an equilibrium 
between ILM creation and destruction is achieved for each energy. This is a very 
low population with no thermodynamical effects but with larger mean energy than 
phonons. This energy is also localized and can be delivered more effectively to the 
bonds that are to be broken. Another example for germanium consists of a series of 
experiments where it was found that Ar plasma ions with energies of 2-8 eV were 
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Fig. 9 Some amplification factors supposing that a potential barrier decreases by AE during half 
the time and increases by the same amount for the other half. The amplification function is simply 
/ = cos\\(AE/ksT) ~ ^tx'p^AE/ksT). A similar process will occur when the barrier oscillates 
while interacting with a moving or stationary breather 


able to anneal defects like the E-center at least two jim below the surface III. On 
the other hand EBD was found to create defects up to a depth of one jim ||T31 . 


6 Conclusions 

It was established that during EBD energetic particles, the product of elastic col¬ 
lisions between 10 keV electrons and residual gas atoms in the vacuum, were the 
primary cause of defects introduced in Ge and Si. High energy electrons interacting 
with the semiconductor directly were found to transfer far less energy, per collision, 
than if the energy transfer occurred through an intermediary atom or molecule. The 
maximum energy transferred via this two-step process was calculated to be approx¬ 
imately 1.3 eV for particles with an atomic mass from 1 to 4 and then diminished 
for heavier particles. This amount of energy, when transferred to a Ge lattice atom, 
is incapable of creating a Erenkel pair but may be sufficient to modify an existing 
defect structure that was previously invisible to DLTS. This conclusion can also be 
drawn if n-Si is used Q. The energies transferred to the germanium lattice by EDB 
is typically of the order of magnitude of intrinsic localized modes. These nonlinear 
localized wave packets have the property of significantly increasing the probability 
of structure changes by temporally lowering the potential barrier for the process. 
Therefore, intrinsically localized modes are very likely to be the cause of the ob¬ 
served phenomenon 








16 


S.M.M. Coelho, J.F.R. Archilla, F.D. Auret and J.M. Nel 


Samples exposed to the conditions of EBD, without deposition (termed EB ex¬ 
posure) did not contain the same defects as the EBD samples except for Eqj? and 
the vacancy-antimony center (V-Sb), Eo.sg. This implies that a necessary condition 
for the introduction of EBD defects was a thin metal layer through which energy 
was transferred to the germanium crystal lattice. The EB exposure defects have not 
yet been identified and may be related to impurities that were accelerated into the 
germanium near-surface region before diffusing deeper into the material, although 
this cannot explain the low defect concentrations observed, especially if the sample 
temperature was allowed to increase during treatment. 
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